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ABSTRACT 

Using VLT/SINFONI, we have obtained repeated AO-assisted, NIR spectroscopy of 
the three central WN6ha stars in the core of the very young 1 Myr), massive and 
dense Galactic cluster NGC3603. One of these stars, NGC3603-A1, is a known 3.77- 
day, double-eclipsing binary, while another one, NGC3603-C, is one of the brightest 
X-ray sources among all known Galactic WR stars, which usually is a strong indication 
for binarity. Our study reveals that star C is indeed an 8.9-day binary, although only 
the WN6ha component is visible in our spectra; therefore we temporarily classify star 
C as an SBl system. Al, on the other hand, is found to consist of two emission-line 
stars of similar, but not necessarily of identical spectral type, which can be followed 
over most the orbit. Using radial velocities for both components and the previously 
known inclination angle of the system, we are able to derive absolute masses for both 
stars in Al. We find Mi = (116 ± 31)M0 for the primary and M2 = (89 ± 16)Mq for 
the secondary component of Al. While uncertainties are large, Al is intrinsically half 
a magnitude brighter than WR20a, the current record holder with 83 and 82 Mq, 
respectively; therefore, it is likely that the primary in Al is indeed the most massive 
star weighed so far. 
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1 INTRODUCTION 

While models maintain that in the early Universe, the first- 
generation of stars were very massive and reached masses be- 
tween 100 and 1000 M© (e.g. iNakamura fc Umemurall200ll : 
ISchaereJbOOSl 'l. it is generally accepted that under present- 
day conditions, relatively fewer massive stars are formed, 
i.e. that the initial-mass function (IMF) is much steeper 
and, more importantly, has a cut-off occu rring around 150 
Mq (IWeidner fc Kroupalliool : lFigejl2005l) . So far, however, 
whenever Keplerian orbits of binary systems are used to 
weigh stars - the only way to obtain reliable, least model- 
dependent masses -, measured masses fall short by almost a 
factor of two with respect to the putative cut-off. Currently, 
stars with the highest known masses are both WN6ha com- 
ponents of the Gal actic WR binary WR20a, with 83 and 8 2 
Mq, respectively l|Rauw et al.ll2004l : ISonanos et al] 12004 ). 
and the 03f/WN6 star in the Galactic binary WR21a with a 
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mmnnum mass of 87 M0 iGamen et al.ll2008l ). Significantly 
more massive stars, however, have so far remained elusive. 

A most remarkable result that has emerged from the 
search for very massive stars is, however, that the high- 
est Keplerian masses are not found among absorption- line 
O-type stars - rat her, masses of these sta r s remain be- 



low ^ 60 M© (e.g. iLamontagne et aL I ll996l : iMassev et al.l 
'20021) but among Wolf-Rayet (WR) stars, more pre- 
cisely among the so-called WN5-7ha (or WN5-7h) stars, 
an extremely luminous and hydrogen-rich subtype of the 
nitrogen-sequence WR stars. Contrary to classical WR 
stars, which are identified with bare, helium-burning cores 
of evolved massive stars, theoretical work confirms that 
these luminous WN5-7ha stars are hydrogen-burning, un- 
evolyed objects l|de Koter et al.l Il997l : ICrowther fc DessartI 
Il998l ) which mimic the spectral appearance of WR stars 
because their high luminosi ties drive dense and fast winds 
ICraefener fc Hamannll2008l ). 

If it is true that the mo st massive stars can be found in 
the most massive clusters l|Batel l2002l : IWeidner fc Kroupal 
l200d ). then it can be expected that the most massive 
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WN5-7ha stars are hosted by the most massive among 
the youngest, least evolved clusters known. NGC3603 is a 
Galactic example of such a cluster, and virtually a clone 
of its more famous LMC counterpart, the supermassive 
cluster R136, at th e core of the giant HIT region 30 Dor 
ijMoffat et al.lll994l l. NGC3603's very core, itself denoted 
HP 97950, contains three extremely luminous WN6ha stars 
l|Drissen et al.lll995h. with stellar luminosities well in ex- 



cess of 10" L, 



__ M (Ide Koter et"all Il997l : ICrowther fc DessartI 

1 19981 '). iMoffat fc Niemelal l|l98i ) found from radial- velocity 
variations in unresolved spectra of HD 97950 that one of 
these WN6ha stars must be a binary with a perio d of 3.772 
days, whic h was confirrned by iMoffat e\~^ l(l985h . In a re- 
cent studv. lMoffat et all (|2004l ) used HST-NICMOS J-band 
photometry to confirm that Al is indeed a double-eclipsing 
binary with that period, while the other WN6ha stars B and 
C did not vary above the ~ 0.05 mag noise level. However, 
star C shows an ex tremely large X-ra y luminosity (Lx = 
several 10 ergs' -1; iMoffat et al.]|2002l ). which is a strong 
indication of binarity given that colliding winds in binary 
syste ms generate copious amounts of hard X-rays (e.g. lUsovl 
1 19921 ). 

The double-eclipsing nature of Al offers the rare oppor- 
tunity to directly measure the mass of an extremely lumi- 
nous binary system by simple, least model-dependent Kep- 
lerian orbits. Here we report the results of a spectroscopic 
monitoring campaign of the three central WN6ha stars in 
NGC3603. This Letter is organized as follows: In Section 
2, we briefly describe the observations and data reduction. 
In Section 3, data analysis is described and results are pre- 
sented. Section 4 summarizes the Letter. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were carried out in service mode at the VLT- 
UT4 under Program-ID P75-0576.D. We obtained repeated 
K-band (1.95 to 2.45 /jm) spectroscopy using SINFONI 
l|Eisenhauer et al.ll2003l : iBonnet et al.ll2004^ with adaptive- 
optics (AO) correction to obtain the highest possible spatial 
and spectral resolution. The field of view was 0.8" x 0.8" 
with a "spaxel" scale of 12.5 mas x 25 mas. Each of our 
three target WN6ha stars was observed individually. 

Total exposure times were 84s per star and per visit, 
each organized in 4 detector integration times (DITs) of 21s. 
Given the brightness of our targets {K ~ 7-8mag) and the 
AO deployment, no dedicated sky frames were taken. Other 
calibrations (dark and flatfield frames, and the telluric stan- 
dard star) were provided by the ESO baseline calibration. 

For mos t of the data reduc tion steps, ESO's pipeline 
was used fcf. lAbuter et al]|2006l ). Standard reduction steps 
were taken. The two-dimensional spectra produced by each 
illuminated slitlet were individually extracted using IRAF, 
and combined into one wavelength-calibrated spectrum per 
star and visit. A main-sequence B-type star was used for 
telluric corrections. To remove the B star's Br7 absorp- 
tion which coincides with the Br7/Hel A2.166 /im emission 
blend of the target WN6ha stars, a Lorentzian was fitted to 
the absorption line and subtracted from the B star's spec- 
trum. Residuals were very small, and proved to be harm- 
less in the subsequent analysis. Finally, science spectra were 
rectified by fitting a low-order spline function to the stel- 
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Figure 1. Montage of average spectra of our three target stars, 
limited to the useful spectral region. The two strongest emission 
lines, Br7/Hel A2.166 /^m and Hcll A2.188 lira are indicated. For 
clarity, the upper two spectra have been shifted by 0.25 and 0.5 
flux units, respectively. 
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Figure 2. Average of two spectra of Al, taken almost at quadra- 
ture ((/) ~ 0.74). Both Br7/Hel A2.166 /^m and Hell A2.188 fiva 
are clearly separated; the (stronger-lined) primary is indicated 
by upward arrows, while the (slightly weaker-lined) secondary is 
indicated by downward arrows. 



lar continuum. The final uniform stepwidth of the spectra 
was 2.45 A/pixel, resulting in a conservative three-pixel re- 
solving power of ~3000, and a velocity dispersion of ~33 
kms^^ /pixel. 



3 DATA ANALYSIS AND RESULTS 

Averaged spectra of the three WN6ha stars are shown in 
Figure [T] For star B, radial velocities (RVs) were measured 
by cross-correlation of the region 21340 to 22140 A, com- 
prising the two emission lines Br7/Hel A2.166/im and Hen 
A2.188/im. Star B proved to show constant RVs over the ob- 
served timescales, with a scatter of ctrv ~ 20 kms~^, which 
was adopted as a-posteriori error (see Figure |4]). Fitting a 
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Figure 3. Power spectrum (upper panel) and spectral window 
(lower panel) for star C. Frequencies were searched from 0.1 to 
1 d— 1, corresponding to periods between 1 and 10 days. The 
strongest frequency peak at v ~ 0.112 d^^ corresponds to the 
adopted period of 8.89 days. 

single Gaussian to the Hell A2.188/^m emission line yielded 
a slightly larger scatter, a — 2Q kms~^, but was used to ob- 
tain the systemic velocity of the star, 7 — (167 ± 6) kms~^. 
Given is the error of the mean, i.e. a/^/{N), where iV = 18 
is the number of data points. 

Al was found to consist of two e mission-line stars, very 
similar to WR20a l|Rauw et al.ll200^ ). An average spectrum 
around quadrature {<f) ~ 0.74) is shown in Figure [S] While 
both binary components show Br7/Hel A2.166/im and Hell 
A2.188/im in emission, the two components do not show the 
same line strengths. In the following, we will hence refer to 
the stronger- lined star as "primary" , and to the weaker-lined 
star as "secondary". 

RVs for star C and for both components of Al were 
obtained by independently fitting Gaussians to the respec- 
tive Hell A2.188/im emission peaks of the two components, 
because this line is most likely less affected by wind-wind col- 
lision in the binary than Br7/Hel. To take into account that 
in the case of Al, the Hell emission peak of the secondary 
moves up and down the flank of the stronger emission of the 
primary, a sloping, linear continuum was simultaneously fit- 
ted. Unfortunately, fitting the secondary's emission peaks in 
Al is only possible at phases were the respective Hell lines 
are well separated, and particularly difficult when the weak 
Hell line of the secondary blends with the strong Br7/Hel 
line of the primary, i.e. during most of the "blue" half-wave 
of the companion's motion. Moreover, due to blending with 
the steep slope of the Br7/Hel emission, the profile of the 
secondary's Hell emission line gets significantly distorted, 
thereby further increasing the error of the fit. Both phe- 
nomena strongly affect the determination of the secondary's 
orbit and hence the mass determination of the two stars (see 
below). RVs for both components of Al are listed in Table[TJ 
and RVs for the primary component of C are listed in Table 

m 

For star C, RVs were subjected to a period search 
base d on the Clean alg orithm and Scargle periodograms 
(see ISchnurr et al] I2OO8I for details). The search yielded a 
value of P = 8.89 ± 0.01 days (see Figure [3}. 

Orbital s olutions for Al and C were fitted to the RVs 
using Elem (jMarchenko et al.l Il994l ) and applying 1 /cr^v 



Table 1. Journal of RVs for both components in Al as obtained 
by Gaussian fits to the Hell A2.188/xm emission lines. The second 
column of each table provides the orbital phase corresponding to 
a period of P = 3.7724 days (from Moffat ct al. 2004) and, as zero 
point in phase, the time of inferior conjunction Eq = 2, 453, 765.75 
(this work). Note that it was not always possible to measure the 
RVs for the secondary; see text for more details. 
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weights to t h e RV s. For star Al, the period reported by 
iMoiffat et al] (|2004l ). P = 3.7724 d, was adopted, and for 
both Al and C the orbital periods were fixed. In Al, RVs 
of both the primary and the secondary were fitted simulta- 
neously, forcing a circular orbit and the systemic velocity to 
be the same for both components. While emission lines in 
WR stars are well known to display some degree of redshift 
wit h respect to the tru e systemic velocity of the star (see 
e.g. ISchnurr et al.|[2008h . both WR components in Al will 
display the same amount of redshift. Our approach is fur- 
ther justified by the fact that all three stars Al, B, and C 
show comparable systemic velocities. 

Due to the low quality of the secondary's RV determina- 
tion during the "blue" half-wave of its orbital motion, the fit 
of the secondary's orbit is dominated by the "red" half-wave. 
This leads to a large RV amplitude for the secondary and, by 
consequence, a large mass ratio. A significantly smaller RV 
amplitude would have been obtained if the systemic velocity 
of the secondary were left as a free fitting parameter. How- 
ever, given the problems with the 'blue half-wave described 
above, we believe that it is more sensible to rely on the red 
half-wave with the smaller error bars and fix the systemic 
velocity. 

Orbital fits for Al and C are shown in Figure |4l the 
respective orbital parameters are given in Tables [3] and U) 
Combining the orbital parameters of b oth componen t s of A l 
with the inclination angle obtained bv lMoffat et al.l (|2004h . 
i — 71°, and using the usual relations, we obtain abso- 
lute masses of Mi = (116 ± 31) Mq for the primary and 
M2 — (89 ± 16) Mq for the secondary component of Al, 
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Table 2. As for Table [T] but for star C, using the period of P = 
8.89 days and, as zero point for the phase, the time of periastron 
passage To = 2,453, 547.11 (this work). 
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Table 3. Orbital parameters for both the primary and the sec- 
ondary component of Al from the combined, weighted fit, forcing 
a circular solution. The inclination angle and the period have been 
adopted from Moffat et al. (2004). 



Parameter Primary Secondary 

P [days] 3.7724 
i [°] 71 
e 

[2,450,000.5+] 3765.25 +0.03 

7 [kms-l] 153 + 12 

K [kms-i] 330 + 20 433 + 53 

o-o-c [kms-1] 42 82 

M [Mq] 116 +31 89 + 16 



respectively. Due to the problems with the secondary's RVs 
described above, the uncertainties on the orbits, expressed 
by the large CTo-c, a nd hence on the mas ses are uncomfort- 
ably large. However, iDrissen et al.l l| 19951 ) have reported an 
absolute magnitude of My ~ —7.5 for Al, i.e. bright er by 
~ 0.5 mag than WR20a, for which iRauw et al.] (|2007l ) find 
Mv = 7.04 ± 0.25. Since Al and WR20a have very sim- 
ilar spectral types and, thus, bolometric corrections, Al's 
brighter magnitude directly translates into a larger bolo- 
metric luminosity, qualitatively consistent with our result 
that Al is more massive than WR20a. However, better data 
are currently being taken to verify the results. 

In star C, the secondary is not visible in the spectrum; 



Table 4. Orbital parameters for the WN6ha (primary) compo- 
nent in star 0. 



Parameter Primary 

P [days] 8.89 + 0.01 

e 0.30+0.04 

[°] 281 + 7 
To [2,450,000.5+] 3,546.61 +0.18 

7 [kms-l] 186 +6 

K [kms-i] 200 + 23 

o-Q-c [kms~l] 27 



hence, we temporarily classify the system as SBl. While the 
intrin sic brightness difference b etween Al and C is ~ 0.5 
mag (jCrowther fc Dessartlll99^ . it is well possible that the 
WN6ha component in C is as massive as the primary in Al, 
where the luminosities of the primary and secondary are 
more similar; for the same reason, a large mass ratio can be 
expected between the primary and secondary component in 
C. 

Star B, on the other hand, di splays neither sign ificant 
RV nor photometric variations (cf. iMoffat et aLll20M ) nor a 
strong X-ray flux. The latter is a strong indication against 
B's possible long-period binarit y. Therefore, B most l i kely is 
a truly single star. Remarkably, ICrowther fc Dessart' ("19981) 
report a luminosity for B which is only very slightly lower 
than that of the (unresolved) system Al. This renders B 
significantly more massive than the primaries in both Al 
and C, and puts it close to (or even in excess of) the pu- 
tative IMF cut-off mass, unless B turns out to be an un- 
resolved, long-period binary as well. Clearly, line-blanketed 
atmosphere models will provide valuable insights into this 
problem once the known, very massive binary systems such 
as Al and WR20a have been used as yardsticks. 



4 SUMMARY AND CONCLUSION 

We have obtained repeated, spatially resolved, AO assisted, 
near-IR spectroscopy of the three central WN6ha stars in 
HD 97950, the core of the young, unevolved and very mas- 
sive Galactic cluster NGC3603. One of the stars, Al, is a 
previously known, d ouble-eclipsing binary w i th an orbital 
period of 3.77 days |Moffat fc Niemelalll984l : iMoffat et al.l 
ll985l : lMoffa~ et al.ll2004l ). while a second star, C, was newly 
identified as a binary in the present study. The third star, B, 
showed constant RVs over the observed time interval, and 
therefore is most likely no t a binary, which is in line with its 
normal X-ray luminosity l| Moffat et al.|[200^ ). 

While in star C only the primary (WN6ha) component 
is visible - the system is therefore classified as an SBl bi- 
nary -, Al consists of two emission-line stars, most likely 
of similar, but hot identical spectral types. From the radial- 
velocity curves of the two com ponents and the k nown in- 
clination angle of the system l|Moffat et ahl |2004 ) , we de- 
rived component masses of Mi = (116 ± 31)Mq for the 
primary and M2 = (89 ± 16)Mq for the secondary, respec- 
tively. Despite the large uncertainties, we consider the pri- 
mary WN6ha component of Al to be the most massive star 
ever directly weighed. 

While the primary component of C might have a mass 
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similar to or even greater than that of Al's primary, it is 
possible that star B, single yet only slightly fainter than 
the combined binary system Al, is indeed the most massive 
member in NGC3603 and, therefore, the most massive main- 
sequence star known in the Galaxy. 
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Figure 4. Orbital solutions obtained from weighted fits for Al 
(upper panel) and C (middle panel), folded in the respective 
phases. In Al, the primary is indicated by filled symbols and 
the solid curve, while the secondary is shown with open symbols 
and the dotted curve. In star C, only the primary is visible. RVs 
of star B arc constant during the observations (lower panel, same 
scale as for star C) with a scatter of cr = 20 kms~^ . The dashed 
line indicates the systemic velocity of the star B, 7 = (167 ± 6) 
kms~^. 



